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Abstract 

We calculate using loop expansion the effect of fluctuations on the structure 
function and magnetization of the vortex lattice and compare it with existing 
MC results. In addition to renormalization of the height of the Bragg peaks 
of the structure function, there appears a characteristic saddle shape "halos" 
around the peaks. The effect of disorder on magnetization is also calculated. 
All the infrared divergencies related to soft shear cancel. 
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I. INTRODUCTION 



Decoration neutron scattering and STM H] clearly demonstrated the Abrikosov flux 
line lattice in low and high T c type II superconductors. There are however important 
differences between the two classes of materials. Ginzburg parameter Gi characterizing 
importance of thermal fluctuations is much larger in high T c superconductors than in the 
low temperature ones. Moreover in the presence of magnetic field the importance of 
fluctuations in high T c superconductors is further enhanced. The lattice melts and becomes 
vortex liquid over large portions of the phase diagram In "strongly fluctuating" 

superconductors, even far below the melting line, corrections to various physical quantities 
like magnetization or specific heat are not negligible. The vortex lattice becomes distorted. 
It is quite straightforward to systematically account for the fluctuations effect on 
magnetization, specific heat or conductivity perturbatively above the mean field transition 
line using Ginzburg - Landau (GL) description 0. However in the interesting region below 
this line it turned out to be extremely difficult to develop a quantitative theory. 
A direct approach to the low temperature fluctuations physics is to start from the mean 
field solution and then take fluctuations around this inhomogeneous solution into account 
perturbatively. Experimentally it is reasonable since, for example, specific heat at low 
temperatures is a smooth function and the fluctuations contribution is quite small. For 
some time this was in disagreement with theoretical expectations. Eilenberger calculated 
spectrum of harmonic excitations of the triangular vortex lattice |§ and noted that the 
gapless mode is softer than the usual Goldstone mode expected as a result of spontaneous 
breaking of translational invariance. The inverse propagator for the "phase" excitations 
behaves as k 2 + const{k 2 x + k 2 ) 2 . It was shown fl9|]T0H that this behavior is directly related 
to the nondispersive nature of the shear modulus cqq and is in agreement with numerous 
experiments. 

The influence of this additional "softness" goes beyond enhancement of the contribution of 
fluctuations at leading order. It apparently leads to disastrous infrared divergencies at 
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higher orders rendering the perturbation theory around the vortex state doubtful. One 
therefore tends to think that nonperturbative effects are so important that such a 
perturbation theory should be abandoned [|TTJ. However it was shown in |L2| that a closer 
look at the diagrams reveals that in fact one encounters actually only logarithmic 
divergencies. This makes the divergencies similar to so called "spurious" divergencies in 
the theory of critical phenomena with broken continuous symmetry and they exactly 
cancel at each order provided we are calculating a symmetric quantity. One can effectively 
use properly modified perturbation theory to quantitatively study various properties of the 
vortex liquid phase. Magnetization calculated using this perturbative approach agrees very 



well with the direct Monte Carlo simulation of |JL3|| . The method was then extended 
beyond the lowest Landau level (LLL) fl4|| . 

In this paper we calculate the effect of fluctuations on the magnetic field distribution, 
structure function of the vortex lattice and compare with existing MC results. 
Fluctuations cause the spread of the peaks in diffraction pattern in a very specific way, 
while height of the peaks is slightly corrected. Effects of fluctuation and disorder on 
magnetization and specific heat are computed. The paper is organized as follows. In 
section II the model and the fluctuation spectrum approximation are briefly reviewed. In 
section III the calculation of the structure function is presented. Section IV contains 
analysis of the result, comparison with MC simulation and some generalizations. In section 
V the distribution of magnetic field is calculated, while effects of weak disorder on 
magnetization and specific heat are treated in section VI. Summary is given in section VII. 



II. MODEL, MEAN FIELD SOLUTION AND THE PERTURBATION THEORY 



Our starting point is the GL free energy: 
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Here A = (—By, 0) describes a nonfluctuating constant magnetic field. For strongly type 
II superconductors (k ~ 100) far from if cl (this is the range of interest in this paper) 



magnetic field is homogeneous to a high degree due to superposition from many vortices. 
For simplicity we assume a = a(l — t), t = T/T C although this dependence can be easily 
modified to better describe the experimental coherence length. 
Throughout most of the paper will use the following units. Unit of length is 



£ = yh 2 / (2m a baT c ) and unit of magnetic field is H C 2, so that dimensionless magnetic field 
is b = B/H c2 . The dimensionless free energy in these units is (the order parameter field is 
rescaled as ip 2 — > ^^i/j 2 ): 
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The dimensionless coefficient is 
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where the Ginzburg number is defined by Gi = |( 32?re 1 2 fp cl 1 ) 2 and 7 = m c /m a b is an 
anisotropy parameter. This coefficient determines the strength of fluctuations, but is 
irrelevant as far as mean field solutions are concerned. 
The second expansion parameter is (see [§,0] for details): 
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If ah is sufficiently small GL equations can be solved perturbatively: 
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It is convenient to represent $cb ^"l? ••• i n the basis of eigenfunctions of operator 7i = 
|(— D 2 — b), TCip n = nb(p n , normalized to unit "Cooper pairs density" 



< \ip n \ >= Jceiid x\(p n \ 2 — = 1, where "cell" is a primitive cell of the vortex lattice. 
Assuming hexagonal lattice symmetry one explicitly has: 
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where = y ^ is the lattice spacing. One finds 
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To order a l h , we expand 

oo 

These coefficients can be found in |14[ . 
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To find an excitation spectrum one expands free energy functional around the solution. 
The fluctuating order parameter field ip is divided into a nonfluctuating (mean field) part 
and a small fluctuation 



ip(x) = + x( x )- 
Field x can be expanded in a basis of quasimomentum eigenfunctions: 
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Instead of complex field xt we will use two "real" fields 0% and A% satisfying 
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where d k = exp[— i9 k /2] where 7& = \j k \ exp[i9 k ]. Within the LLL, the eigenstates are 
A^, 0%, while the eigenvalues (in two dimensions; in three dimensions simply plus are 



e A = a h e\ = a h ^-1 + -f3 k - ^\lk\^j 
to = a h e = a h ^-1 + ^(3 k + jjllkl^ 



(12) 
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where €a, to ar e dependent on two dimensional vector k. Higher order corrections and 



higher Landau levels eigenstates and eigenvalues can be found in |T^]. With spectrum of 
excitations and expansion of solutions of GL equations in ah one can start the calculation 
of correlators to any order in to. 

III. STRUCTURE FUNCTION OF THE VORTEX LATTICE 

In this section the structure function is calculated to order u (harmonic approximation) 
within the LLL, namely neglecting higher corrections. We discuss these corrections in 
the next section. First we calculate the density correlator defined by 

3(z, *a) = (p( x , ^s)p(x + z, x 3 + z 3 )) = (p(x)p(y)) x (13) 
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where <> x indicates average over x (which means here over the unit cell) and p = 
The correlator is calculated using Wick expansion: 

S = S rnf + uSfl uct . (14) 

The first term is the mean field part, while the second term is the fluctuation part. The 
mean field part is simply 



S mf = (Mx)\ 2 My)\ 2 ) x . (15) 



The structure function is the Fourier transform >S(q, 0) = / dze l<i ' z S(z, z% = 0). Within the 
LLL, $(a;) = (jfc) ^ (fix) an d the mean field part of the structure function becomes: 
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where we made use of formulas ..and function <5 n (q) defined in Appendix. This is just sum 
of delta functions of various heights at reciprocal lattice points. 
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The fluctuation part contains four terms (diagrams) Si, ...S4. The first term is 



Si(z,z 3 ) = i _L_ ^$(x)$(y) f ki <r k d!<p?{x)<p?{y) 
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< OfeO™ > and < A^Af > are propagators: 
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To calculate structure functions we will need only the z 3 = correlator: 



Si(z,0) = 
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Within the LLL approximation it simplified to: 
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The first fluctuation correction term to structure function can be evaluated as follows: 
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where formulas of Appendix were used. Q is the integer part of q, k is the fractional part 
of q: q = k + nidi + n 2 d 2 = k+ Q (see Appendix for the definitions of d\, d 2 ). The 
second fluctuation correction term is 
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S 2 (z, z 3 = 0) is equal to (in the LLL approximation): 
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The third term is 
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and within the LLL at z 3 = is equal to: 

s 3 (z, o) = ^J-p (|$(x)| 2 jf ^(y)</W 

+ (x < — ► y). 
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Consequently the correction to the structure function is: 
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The final term is from the vacuum renormalization contribution. The shift v in 
ip(x) = v<j>(x) + x( x ) is renormalized, that is, to one loop order, v 2 = v 2 + uvf, where 
v o = fj- One can find v \ by minimizing the effective one loop free energy: 
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e A (k) = -a h + 2v 2 p k - v 2 ^^ 
e (k) = -a h + 2v 2 p k + v 2 \-f k \. 



Straightforward calculation gives: 
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The last contribution to the one loop correction to the correlator is therefore: 
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The sum of all the four terms can be cast in the following form: 
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= -28.52756. 



The reason we write the sum in such a form will be explained in the next section. 
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IV. ANALYSIS OF THE RESULT, COMPARISON WITH MC SIMULATIONS 



Although all of the four terms Si, S 2 , S 3 and £4 are divergent as any of the peaks is 
approached, k — > 0, the sums Si, S2 and 53, S4 are not. We start with the first two: 
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5 1 (q,0) + 5 2 (q,0) = ^exp 
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where /i(q) defined in eq.(^) contains a function \{k x k y + k x Q) + 6 k . When k — ► it 



can be shown that + 6 k = O (k 2 ' 2 ) , thus f (A;^ + k x Q) + 9 k -> k x Q , and 



1 _ CQs( WkxQ + 



(k x Q) 2 . Hence it will cancel the 1/A; 2 singularity coming from 



p-pjy- Thus /i(q) approaches const. + const. ■ when Q 0, and approaches 

const. + const. ■ k 6 when Q = 0. 

Similarly the sum of S^q, 0) and S^q, 0) is not divergent, although separately they are. 
Their sum is. 
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Now we compare our results with numerical simulation of the LLL system in ||13|| . The 
general shape of the structure function in the vicinity of a peak (see Fig. lb) and the data 
near the origin according to a MC simulation of the same system within the same LLL 
approximation in ref. [13[] (Fig. la) are qualitatively same pattern. It is easier to compare 
using rescaled quasimomenta: q — > qVb, k — > ky^. We get 
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where /i(q), / 2 (Q) and / 3 are defined in eq. (|32"D , but with 6=1 (for example, 
fs = —28.5275) and the region of integration in the formula rescaled to the cell with 
di,d2 being the reciprocal lattice basis vectors 

2vr / 
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Furthermore we define s(q) which is used also in [iU 

Qf n\ ( ah Y 471-2 [ q2 l ^ 
5(q,0)=^-j — exp^ T j S (q) 

S (q) = 5 n (q) + ^M^) [/l(q) + 5n (q)f 2 (Q) + 5 n (q)f 3 ] 
For reciprocal lattice vectors close to origin the value of /^(Q) are: 
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Table 1. 

Values of /2(Q) with small ri\n%. 



n\n 2 


0,1 


1,0 


1,1 


1,2 


0,2 


2,2 


1,3 


MQ) 


-46.19 


-46.19 


-46.19 


-63.165 


-64.35 


-64.35 


-73.825 



In ref. |13[], the material parameters describe YBCO: T c = 93 K, 

dH c2 {T)/dT = -1.8 x 10 4 Oe/K, 7 = 5, and k = 52. At T = 82.8K, H = 50 kOe. This 
leads to following dimensionless parameters Gi = 2.04 x 10~ 5 , uj = 0.056, = 0.039904. 



However as discussed in |1J] effective expansion parameters are |^ = 0. 222 68 and 



c ~i- 3 / 2 _ 

2^/2^ — = 0.79328, both less than one. |£ is the parameter for the expansion of the 



classical solution. The factor 6 comes from a fact that due to hexagonal, thus only 6th, 
12th, etc. Landau levels appears in perturbation expansion. ^^J^J * s the parameter for 
the fluctuations (loop) expansion and is just barely less than one here. It justifies the 
quantum correction of the formula using perturbation expansion. The numerical factor in 
front of the fluctuation correction in this case is 

(3 A buj {a h y i/2 
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In a finite size sample, 8 n (q) is equal to j^jz when q lies on the reciprocal lattice 
midi + 7B2CI2, otherwise it is zero. Because L x L y = N^2tt (N^ is the number of vortices of 
order 100 only in the MC simulation), thus it is equal to The normalized structure 
function (s n (0) = 1, as it was used in |13|]) is: 
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The correction to the height of the peak at Q, tj^^t/2(Q), is quite small. We find the 



height of peak away from origin found in the MC simulation |L3[ are typically smaller than 
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ours, while around the peaks are larger than analytical. It may be due to finite size effect 
or finite samplings of MC calculation. In the MC calculation part of the peak might 
"belong" to a neighboring pixel. We plot the correction to the non-peak region on Fig. lb 
and find that the theoretical prediction has roughly the same characteristic saddle shape 



"halos" around the peaks as in ref. |13| , Fig. la on which all the peaks were removed (so it 
is different from Fig.2a in fl3| on which only the central peak was removed). 
We can extend our formula to higher orders which will include also the HLLs (higher 
Landau levels). To next order of ah, we should include $1 in $, $ = (a^) 1 / 2 [$o + ct/i^i], 
consider eo(k),€A(k) to next order a\, and €£,(£;), e^(/c) to order a^. It is straightforward to 
do it. 



V. FLUCTUATION OF MAGNETIC FIELD 



Another quantity which can be measured is the magnetic field distribution. In addition to 
constant magnetic field background there are 1/k magnetization corrections due to field 
produced by supercurrent. To leading order in it is given by m{x) oc <p ^ > (for example, 
see ref. [0). (p(x)) can be calculated using the following equation, 



(40) 
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= |$(x)| 2 + ($*(a:)x(a:)) + (${x)x*(x)) + {x(x)x*(x)) 

= ^\Hx)\ 2 + (x(x) X *(x)). 
Pa 

Using eq.(prS[) and eq. (]TT|), and considering only x 3 = 0, one obtains: 
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However, as pointed out in Sec. Ill the coefficient v in ip(x) = v<p(x) + x( x ) is 
renormalized to one loop order, v 2 = v 2 Q + uw 2 , with v\ given in eq.(^). Thus we need to 
add a term, iov\ [0(x)] 2 to eq.([4C)|). 
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Its Fourier transform p(q) = J dze iq ' z (p(x, 0)) can be easily calculated: 
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q . n*% , m f 2 



I— H , 

l /3^ 16vr 2 A 

/ ( ik x q N 

+ exp 



cxp 



ik x q N 



'2/3* + |7fc| 







eo(k) 



2/3 fc - f-yjfc J 



, b j \ p ;;v^(k) j 

Performing integrals and rescaling the quasimomenta again, one obtains: 
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The function /^(Q) appeared in eq.( 



VI. DISORDER EFFECT ON MAGNETIZATION AND SPECIFIC HEAT 



One can introduce weak disorder by adding a quadratic term in eq. (@) ||, 



Af = / d 3 xa(x)\^ 



(44) 

Loosely speaking it represents local variation of temperature. For pointlike defects one can 
assume that the correlation of a(x) is << a(x)a(y) »= W5(x — y), << a(x) »= 0. 
Before the disorder average we calculate the free energy — T In Z with 



Z = J VcVcox\) {-- /[V>*, V] - / d\niU)\t f 



(45) 



If W is very small, we can calculate Z by perturbation theory in W. To the second order 
Z is given as 

If,, . s 1 



Z = Zz 
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where Zq is the free energy without disorder and it had been obtained in ref. [ 14| . Thus 
the free energy with disorder is 

F = -T In Z = F + AF (47) 

r q/(x) T ft 

= F + T < p(x) > - — / / [< p(x)p(y) > - < p(x) >< p(y) >} a{x)a{y), 

Jx LU ZUJ Z Jx Jy 

where F = —TlnZ . Averaging free energy over disorder one obtains: 

TW r 

F = F - -^jj x [< P(x)p(x) >-< p(x) >< p(x) >} (48) 

TWV ~ 
— Fo ^-p-u>Sfi uc t(0). 

From eq. (p2[) , one finds that 5/^(0) = -0.18619^ . Hence the energy density 

1 /2 

difference due to disorder is T = ^f- = — 0.0931 TW Jp h A ~ ' • Since u = \j2Gm 2 t T = ca]/ 2 b 
with c = — °f 9 J i , 1T ^ ■ The disorder effect on magnetization and specific heat are 

Am =- d -w = - c - r" 12 ) < 49) 

respectively. 

VII. CONCLUSIONS 

To conclude, we have calculated the effect of fluctuations on the structure function of the 
vortex lattice and compared it to existing MC results. In addition to renormalization of 
the height of the Bragg peaks, there appears a characteristic saddle shape "halos" around 
the peaks as found in ref. ||13|| . The calculated fluctuations contribution to the magnetic 



field can be more easily observed in low temperature strongly type II superconductors. 
Finally the predicted dependence of magnetization and specific heat on disorder via 
fluctuations also can be experimentally studied. 

Correlations in flux lattices can be experimentally measured using neutron scattering as 
well as some other more exotic methods like muon spin relaxation, electron tomography, 
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scanning SQUID microscopy etc. [p]-|3|, P^6|JT7|| ■ It would be interesting to detect the effect of 
fluctuations given in the present paper directly from experiments by subtracting the 
"background" of the well known mean field correlator. The calculations show that infrared 
divergencies naively expected in all of the physical quantities calculated above due to 
"supersoft" shear modes in the large k limit cancel. This strengthens the view that the 
loop expansion is a reliable theoretical tool to study the fluctuations effects in vortex 
lattice below the melting point. 
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APPENDIX A: 



In this appendix, we present some basic formulas used in the calculations. The basic 
matrix element is: 
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2^ 



/ dxip(x)ipl(x) exp[-ix • q] = A qik exp 

J cell 



7ri . 2 x q 2 iQxQy , ik x q y 



(Al) 



The Kronecker delta is defined by: 



A q , k = A(q - k) 



1, if q = k + nidi + n 2 d 2 



(A2) 



0, otherwise 

where integers n x = ^di • (q — k) and n 2 = ^rd 2 ■ (q — k). Here d\, d 2 are the reciprocal 
lattice basis vectors 



(A3) 



which are dual to di = (a/y/b, 0), d 2 = (^^, an d a = \Zt§- Integrating over the 
sample area A , one obtains: 



/ dxip^ip^x) exp[— ix ■ q] = 47r 2 <5„(q — k) exp 

J A 



m f 2 \ 



x exp 



'46 26 6 



(A4) 



where <5 n (q — k) is defined as 5 n (q — k) = J2 mi m 2 <Kq _ k — midi — m 2 d 2 ) 
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Figure captions 



Fig. la 



Structure factor from the MC simulation of ref. [13]. The peaks at reciprocal lattice points 
are removed. 

Fig. lb 

Fluctuation correction to structure factor of the Abrikosov vortex lattice, eq. fl38l) . The 
peaks at reciprocal lattice points are removed. 
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